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Abstract. More than a decade ago it was noticed that 
an unexpectedly large value of the measured cosmic an- 
tiproton flux at low kinetic energies could be interpreted 
in terms of a neutralino- induced component. An overpro- 
duction of low energy antiprotons seems however to be dis- 
favoured by recent data from the Bess experiment. The 
strategy we propose here is to focus instead on the high 
energy antiproton spectrum. We find cases in which the 
signal from neutralino annihilations in a clumpy halo sce- 
nario can be unambiguously distinguished from the cosmic 
ray induced component. Such signatures can be tested in 
the near future by upcoming experiments. We investigate 
also the possibility that antiprotons have a finite lifetime. 
We find cases in which the neutralino-induced flux is con- 
sistent with present data for Tp as low as O.f 5 Myr, in clear 
violation of a recent bound derived from the antiproton 
cosmic ray flux. 



1. Introduction 

The idea of exploiting cosmic antiprotons measurements 
to probe unconventional particle physics and astrophysics 
scenarios is certainly not new. The investigation on ex- 
otic antiproton sources was stimulated by early reports 
(Golden et al. 1979; Buffington et al. 1981) of unexpect- 
edly large values for the antiproton to proton ratio in cos- 
mic rays. 

A possibility which has been discussed in some details 
is that a significant contribution to the antiproton flux 
may be given by pair annihilation of relic neutralinos. The 
lightest neutralino, plausible the lightest supersymmetric 
particle in the Minimal Supersymmetric extension to the 
Standard Model (MSSM), is one of the leading dark mat- 
ter candidates. Its coupling with ordinary particles has 
the strength of weak interaction. This is exactly what it 
is needed to provide a relic density of the right order for 
a flat Universe. It implies as well that, if one invokes a 
population of relic neutralinos in the halo of galaxies to 
solve the dark matter problem, these neutralinos can anni- 
hilate in pairs into standard model particles, which might 
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eventually be detected. In particular it has been shown 
that the neutralino-induced antiproton flux can be at the 
level of the background flux from cosmic rays (see e.g. 
Silk & Srednicki (1984); Stecker et al. (1985); Bottino et 
al. (1998); Bergstrom et al. (1999b)). The question to ad- 
dress is of course if there are distinctive features in this 
exotic component so that it can be singled out unambigu- 
ously. 

The ordinary cosmic ray induced antiprotons are gen- 
erated in collisions of primary particles with the inter- 
stellar medium. The main source are pp collisions which, 
for kinematical reasons, produce a characteristic spectrum 
peaked at a kinetic energy of about 2 GeV and rapidly de- 
creasing at low energies. The component from neutralino 
annihilations may not fall as fast; in some cases its max- 
imum is actually in the low energy region. This is the 
signature which was proposed more than a decade ago 
in connection with early antiproton measurements (Silk 
& Srednicki 1984; Stecker et al. 1985). The more recent 
(and probably more accurate) data from the Bess exper- 
iment (Matsunaga et al. 1998; Orito 1998) indicate un- 
fortunately that the low energy flux may not be so high 
as previously thought. It actually seems to be at a level 
that is in very good agreement with a standard prediction 
for the secondary flux if one takes into account pHe colli- 
sions and energy losses during antiproton propagation (see 
Bergstrom et al. 1999b, hereafter Paper I). If future mea- 
surements will confirm Bess data, this will lead inevitably 
to an impasse in this neutralino detection technique: if 
the measured flux is consistent with some estimate of the 
background, there is no way in improving the signal to 
background ratio. 

There are however signatures of a neutralino-induced 
flux which are alternative to the one we have just dis- 
cussed. The strategy we propose here is to search for an 
exotic component in the high energy antiproton spectrum, 
above the expected maximum of the background. This 
kind of investigation is motivated by the fact that the sec- 
ondary antiproton flux at high energies is predicted with 
great confidence. It is remarkable that in this energy re- 
gion nearly all estimates in the literature are consistent 
with each other. The three main ingredients to compute 
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this part of the antiproton spectrum, unUke in the low 
energy range, are actually very well established: data on 
scattering of protons on nuclei are sufficiently aboundant, 
the spectral index of the primary proton flux is known 
with good accuracy and, finally, cosmic ray measurements 
fix quite strictly the 0.6 power law scaling of the diffu- 
sion coefficient with rigidity. It follows then that at least 
the shape, if not the normalization, of the interstellar an- 
tiproton flux at high kinetic energies is signiflcantly con- 
strained. Moreover in this regime solar modulation, which 
introduces a further factor of uncertainty at low energies, 
does not play a main role. 

It was shown in Paper I that in the canonical sce- 
nario in which neutralinos are assumed to form a perfectly 
smooth "gas" of dark matter particles extended through- 
out the halo of the Galaxy, the induced antiproton flux is 
not sufficiently enhanced to cause significant distorsions 
to the high energy secondary spectrum. As we will show 
in the next Section a scenario with a clumpy halo can be 
much more favourable. 

We have singled out three classes of examples in which 
the signal we propose can be unambiguously distinguished 
from the secondary component. In the first two the main 
ingredient is clumpy neutralino dark matter, in the third 
we consider also the exotic possibility that antiprotons 
have a finite lifetime. Before describing these three cases, 
we introduce the formalism needed to give predictions for 
antiproton fluxes from clumps of neutralino dark matter 
in the Milky Way. 



2. Neutralino-induced antiproton flux in a 
clumpy Iialo scenario 

2.1. General discussion 

There are reasons to question whether the distribution 
of dark matter in galactic halos has to be regarded as 
smooth on all scales. Several theoretical models predict 
that dark matter may aggregate in high density substruc- 
tures, "clumps" of dark matter (see e.g. Silk & Szalay 
(1987); Silk & Stebbins (1993); Kolb & Tkachev (1994)). 
From the phenomenological point of view, as the annihila- 
tion rate depends quadratically on the neutralino number 
density locally in space, the presence of such clumps in the 
Galaxy may significantly enhance the chances of detect- 
ing neutralino dark matter. If we postulate that at least a 
fraction of the dark matter in the Milky Way is clustered 
in high density regions, we can associate to a given clump 
of density profile pdird) and located at the position Xci, 
the following antiproton source term: 




j d\,i{p,i{r,i)f 5^{x-xa) ■ (1) 



The first terms on the right hand side of this equation 
are related to the particle physics dark matter candidate. 
In our case, UannV is the total annihilation rate of non 
relativistic neutralinos, while m-^ is the neutralino mass 
(which enters in the equation with a -2 power as the dark 
matter density divided by is the neutralino number 
density). We then sum over the annihilation final states 
F which can give antiprotons in a fragmentation or decay 
process, namely heavy quarks, gluons, gauge and Higgs 
bosons. For each of them, and dN^ /dT are, respec- 
tively, the branching ratio and the fragmentation func- 
tion. The hadronization of all final states has been simu- 
lated with the Lund Monte Carlo program Pythia 6.115 
(Sjostrand 1994). 

o'annV and B^ are computed in the framework of the 
MSSM as described in Paper I. Both the value of the an- 
nihilation cross section and the relative importance of dif- 
ferent final states can vary largely in different regions of 
the MSSM parameter space. The phenomenological ver- 
sion of the MSSM we use is described in Bergstrom & 
Gondolo (1996). We just remind here few definitions we 
need in the discussion in Section ^ The lightest neutralino 
is the lightest mass eigenstate obtained from the super- 
position of the superpartners B and of the neutral 
gauge bosons, and of the neutral CP-even Higgsinos 
and The neutralino is usually defined Higgsino-like if 
it is mainly a linear combination of Hi and i?2 > while it is 
called gaugino-like if it is mainly given by the superposi- 
tion of the B and interaction eigenstates. The formal 
classification is done introducing the gaugino fraction Zg. 

The remaining terms in the expression for the antipro- 
ton source function, Eq. (^ , are related to the neutralino 
density distribution; we have actually assumed that the 
extension of clumps is much smaller than galactic scales. 
The source is hence treated as pointlike: this is not strictly 
necessary but simplifies the discussion which follows. 

The propagation of antiprotons in the Galaxy is 
treated with a two zone diffusion model. Assuming a cylin- 
drical symmetry for the Galaxy, a solution to a transport 
equation of the diffusion type can be derived introducing a 
Bessel- Fourier series to factorize out the radial dependence 
and then solving a one dimensional differential equation 
in the vertical direction. Details on the model and an an- 
alytic solution in case of a pointlike source are given as 
well in Paper I. As the diffusion equation we apply is lin- 
ear in the antiproton number density N , we can consider 
each source separately and then sum the contributions to 
compute N. 

To characterize clumps, it is useful to introduce the 
dimensionless parameter 

<5=i^^!;jii^44^ (2) 

Po J a-^rci Pci[rci) 

which gives the effective contrast between the dark matter 
density inside the clump and the local halo density pq. We 
introduce this definition in Eq. (|l|) and then, switching to 



Piero UUio: Signatures of exotic physics in antiproton cosmic ray measurements 



3 



cylindrical coordinates, we derive the coefficients of the 
Bessel-Fourier series for the source function (see Eq. (18) 
in Paper I). We find: 



'(z) = (cTannW) 



1 



^dN' . _ 1 



X p 



■B" po Ma 6 ■ 



■ 2 T ^ 2. fcx -^fc (''^^W' ) cos(/e6'd)<5(z-Zci) 
Rh Jk+i [vD \ Rh, 



(3) 



where Md is the mass of the dark matter source we 
are considering. This formula has to be substituted into 
Eq. (21) and (23) of Paper I to derive the antiproton num- 
ber density and flux. It is useful to factorize out in the 
expression for the flux the dependence on the MSSM pa- 
rameters. We define, in analogy to the expression for a 
smooth halo scenario, 



*?'(r„,T) = (a. 



F 



Po 5 



dT 



where the coefficient Cpj-^p, which has the dimension of a 
length to the power —2 times a solid angle to the power 
— 1, contains the dependence of the flux on the neutralino 
distribution locally in space and on the parameters in the 
propagation model. There is a one to one correspondence 
between the antiproton flux computed summing the con- 
tributions from the smooth halo scenario as done in Pa- 
per I, and that from a single clump, Eq. ^; one has just 
to replace the coefficient Cprop in Eq. (40) of Paper I with 



Mci 6/po. To facilitate the comparison we will plot 



in units of 10^^ cmsr ^ kpc 



2.2. Antiproton flux from single clumps 

We first suppose that there are few heavy clumps and 
study their individual influence on the result for the 
antiproton flux. Given the factorization introduced in 
Eq. (||), we just analyse the coefficient Cp'^p and its be- 
haviour as a function of the position of the source in 
the Galaxy. Intuitively we expect that sources which are 
nearby give the largest contributions to the flux. In case 
of a propagation model which can be regarded to be as for 
an infinitely large diffusion box, we would find that Cp^^p, 
rather than being a function of the three coordinates Xd , 
depends just on the distance from the observer, the line 
of sight distance L. The boundary conditions in the prop- 
agation model, i.e. free escape at border of the diffusion 
region, break strictly speaking this symmetry; their effect 
is however not very severe for sources which are not too 
close to the border of the diffusion zone (in the vertical 
direction not closer than 0.8-1 kpc). 

In Fig. |l| we plot Cp^^p as a function of the distance L, 
choosing T=l GeV and the set of canonical parameters for 
propagation model introduced in Paper I (diffusion zone 
height hfi ~ 3 kpc, constant term in the diffusion coefh- 
cient = 6T0^^ cm^s~^ and rigidity scale in the diffusion 
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Fig. 1. Cpi-op at T=l GeV and as a function of the dis- 
tance from the observer L. 



coefficient Rq = 3GV). The function we have drawn was 
computed fixing Zd and moving away from our position 
in the Galaxy in the direction of the galactic centre.^] As 
can be seen, we find at large distances a very accurate 
exponential scaling, while a sharp cusp appears below a 
distance of few kiloparsecs: the intuition that if the source 
happens to be close to us it gives a large contribution is 
confirmed. To show how accurate the assumption is that 
Cp[op is essentially just a function of L, we plot in Fig. || 
its isolevel curves in the plane Zd = 1 kpc. We have chosen 
to normalize CpJ^p to 1 at 1 kpc above our position in the 
Galaxy. The isolevel curves are basically circumferences, 
which are just slightly deformed by our choice of fixing 
the antiproton number density to be zero at the border of 
the diffusion zone in the radial direction, i.e. r = 20 kpc 
(thick solid line in the figure). 

We will not examine here the dependence of Cp'^p, 
and hence of the flux, on the parameters which define the 
propagation model: the results are analogous to what it 
was derived for Cprop in the smooth halo case. We just 
mention that introducing a galactic wind fixes a preferred 
direction of propagation and breaks therefore the scaling 
with distance. 

As an example we consider the contribution to the 
antiproton flux given by a clump of neutralino dark 
matter defined by the same choice of parameters as in 
Bengtsson et al. (1990), where dumpiness was introduced 
in connection with the neutralino induced 7-ray signal 



^ Details on how to implement the numerical formula to com- 
pute Cprop given in UUio (1999). 
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Fig. 2. Isolevel curves for Cp^^p in the plane Zd = 1 kpc. 

Cpl^p has been normaUzed to 1 at 1 kpc above our position 
in the Galaxy, ro is our galactocentric distance. 



with continuum energy spectrum. We fix hence Md ~ 

10^ Mq and S ^ lO"'*; the prcfactor MdS is then about 
1.3 • lO^pokpc^, where we assumed as local halo den- 
sity pQ = 0.3GeVcm~^. Comparing the coefficient Cp'^p 
with the analogous quantity in a smooth halo scenario 
(with the same choice of propagation model parameters 
Cprop(T = IGeV) ~ lO^^cmsr-i) we find that the an- 
tiproton flux from this single dark matter clump can be 
at the level (or much higher) than the sum of the contri- 
butions from the whole dark matter halo if this source is 
within about 4.5 kpc. 

We might also consider an opposite approach. It is well 
established that a dark mass of at least 2 ■ IO^Mq is con- 
centrated within 0.015 pc at the galactic centre (Eckart 
& Genzel 1996), forming probably a massive black hole, 
possibly the astrophysical object which is called Sgr A*. 
Such accretion of matter might be associated to a region 
where the density of neutralino dark matter is enhanced 
as well. In an extreme scenario (Berezinsky et al. 1992) 
the potential well of a very steep dark matter halo profllc 
(~ 1/r^'*) is the seed for the formation of the black hole 
itself. It is also intriguing that an excess in the high en- 
ergy 7-ray flux from the galactic centre region, which has 
been found from the analysis of EGRET data (Mayer- 
Hasselwander et al. 1998), can be explained in terms of 
nciitralino annihilations if an appropriate enhancement of 
the neutralino density is present there (UUio 1999). Re- 
gardless of its possible origin, we can estimate how large 
the accretion of neutralinos at the galactic centre should 
be to give a measurable primary antiproton flux. Assum- 



ing that the galactocentric distance is 8.5 kpc, we flnd that 
the flux induced by a source at the galactic centre is pro- 
portional to C^5.op(8.5kpc) ~ 4.5 • lO^^cmsr-^kpc"^ (at 
T=l GcV ). If we require for instance that its contribution 
should be at least one half of the total flux in a smooth 
halo scenario, we find that MdS should be at least of the 
order of W^'^Mq. 

2.3. Collective effects of elumpiness 

There is the possibility that a large fraction of the dark 
matter mass is in clumps, in the extreme case all of it. To 
avoid violating dynamical constraints, clumps should be 
light, with masses probably less than Md ~ lO** — lO^M©. 
If one could deal with a model that gives some accurate 
prediction for the masses of clumps and their distribution 
in the Galaxy, it would be possible to exploit the approach 
of the previous paragraph and estimate the antiproton 
flux by adding the contributions from individual soiirccs. 
As very little is known about the inherently non linear 
problem of generating dark matter clumps, it seems more 
reasonable to follow a probabilistic approach. 

Let / be the fraction of dark matter in clumps and 
Nd the total number of clumps, all roughly of about the 
same mass and overdcnsity. We can dcflnc a probability 
density distribution of the clumps in the Galaxy which in 
the limit of large /, to fulfill dynamical constraints, has 
to follow the mass distribution in the halo. In a Cartesian 
coordinate system with origin at the galactic centre, the 
probability to find a given clump in the volume element 
d^x at position x is: 



Pd{x) d^x = p{x) d^x . 



(5) 



We introduced here A'/^, the total mass of the halo, so that 
Pd has the correct normalization J pd{x)d^x = 1. 

The antiproton source function in the volume element 
d^x at the galactic position x is then: 



g|'(r, X)d'x^ (a^nnV)^ J2 

rriv ■^-^ dl 

X p 

■NdPd{x) j d^Td {Pd{rd)f d^x 

= /(5— ^(aan„^^)2^-T^S d^X, 



(6) 



while the antiproton flux is: 



$;-(ro,T) = (aan„^;)^^B^/^ 



m 



.2 



d^x CpJ.Qp(r, x) 



Po 



dN^ ^^^^ po' 



-Knn.)E^^"/^^Cop(^)- 



(7) 
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In the equation above, Cp^.^p{T) may actually be com- 
puted in the same way as the corresponding coefficient 
for a smooth halo Cprop(7'), just by replacing in the for- 
mer {p{x)/pof' with p{x)/po. The two coefficients have 
similar behaviours; for our canonical diffusion parame- 
ters and an isothermal sphere as dark matter halo profile, 
Cprop ~ 0.75 • Cprop for most of the kinetic energies of 
interest in our problem. The antiproton flux in the many 
small clump scenario can then be obtained by scaling the 
result in the smooth halo case by roughly 0.75 •/ 5. There 
is quite a large freedom in the choice of / S, however one 
has to worry about not violating existing experimental 
bounds (Bergstrom et al. 1999a). In all the results in the 
next Section we have been careful to check that the mod- 
els we propose do not induce an overproduction of both 
diffuse 7-rays and of cosmic positrons. 

3. Signatures in the high energy antiproton 
spectrum 

We have provided two schemes in which the antiproton 
flux can be sensibly enhanced with respect to the results in 
the smooth halo scenario. We propose here three cases in 
which the signal from neutralino dark matter annihilations 
in a clumpy halo can distort the high energy cosmic ray 
flux, giving a very clear signature of the presence of an 
exotic component. 

Two experiments have measured the antiproton flux 
above a kinetic energy of few GeV (Golden et al. 1979; Hof 
et al. 1996) and their data are in contradiction with each 
other. A new set of experiments, the Caprice experiment 
(Boezio et al. 1997) and the space-based Ams (Ahlen et 
al. 1994) and Pamela (Adriani et al. 1995), should give 
much more abundant data in the near future. We will 
check that our predictions are consistent with the data in 
the low energy regime (below 3 GeV) from the Bess 97 
(Orito 1998) and Bess 95 (Matsunaga et al. 1998) flights, 
which have given a first hint on the actual shape of the an- 
tiproton spectrum. The prediction for the background for 
our canonical set of parameters actually provides a very 
good fit of the data (see Paper I; our background predic- 
tion has been confirmed by Bieber et al. (1999) and is con- 
sistent in the high energy region with most the previous 
estimates in the literature). In the first two cases below 
we allow a little room for a neutralino-induced component 
by lowering the normalization of the primary proton flux 
by 1 cr. 

For simplicity we will focus on the many small clumps 
scenario and quote in each case the value of the parame- 
ter f6 we are considering. The sample of supersymmetric 
models on which this analysis is based is the same as in 
Paper I. As in the smooth halo case we restrict to those 
MSSMs for which the neutralino has a cosmologically in- 
teresting relic density, e.g. 0.025 < V.^h'^ < 1. 

To compare with Bess data we have to take into ac- 
count solar modulation. We relate the flux measured at 
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Fig. 3. Distorsion on the secondary antiproton flux (line b) 
induced by a signal from a heavy Higgsino-like neutralino. 

the top of the atmosphere to the interstellar flux applying 
the force-field approximation by Gleeson & Axford (1967) 
with solar modulation parameter (pp ~ 500 MV (value 
suggested in the analysis of the Bess collaboration). 

3.1. Broadening of the maximum in antiproton flux 

The prediction for the background secondary flux has a 
well defined maximum at about 2 GeV. We check if it 
is possible to introduce a distorsion in the spectrum by 
adding a component which widens the shape of the max- 
imum. We will further require that this is not just some 
local distorsion which might be mimicked by introduc- 
ing for instance reacceleration effects (Simon & Heinbach 
1996). In our extreme case we impose that the antiproton 
flux from neutralino annihilations must exceed the esti- 
mate for the background by at least one order of magni- 
tude at very high energies, say about 50 GeV. To obtain a 
wider maximum one has to add a neutralino signal rather 
sharply peaked at an energy higher than the maximum in 
the background; its spectrum should then decrease rapidly 
at low energies. The effect we are searching for might be 
produced by high mass neutralinos with negligible branch- 
ing ratio into bb or tt, which is the case e.g. for a very 
pure heavy Higgsino-like neutralino. In Fig. ^ we show 
two examples compared to our standard background es- 
timate (solid line denoted with the letter b). Model 1 is 
a 964 GeV Higgsino {Zg = 0.0027) with flux rescaled by 
fS~ 4180; it is the model in our MSSM sample which 
gives the highest possible flux at 50 GeV compatibly with 
the background normalization we have chosen (dashed line 
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Fig. 4. Distorsion on the secondary antiproton flux (line 
b) induced by a neutralino signal with flat spectrum. 

in the figure) and with the requirement that added to the 
former it should give a good fit of the Bess data. Model 2 
is instead the case in this class of examples which is associ- 
ated with the lowest rescaling, f 6 = 1180; it is at the same 
time the model with the smallest mass, an Higgsino-like 
neutrahno {Zg = 10"'') with = 777 GeV. A further 
reduction in the background or a less pronounced overpro- 
duction of antiprotons at high energies, allow lower values 
of both and / 6. 

3.2. Break in the energy spectrum 

We focus in this case on rather flat neutralino-induced 
spectra. If one adds such a component to the background, 
the peak and the low energy tail in the antiproton flux 
are not sensibly affected, but a severe break may occur 
in the high energy region. A flat signal spectrum is gen- 
erated overlapping a low contribution from a quark anni- 
hilation channel to a gauge boson contribution at higher 
energies. In Fig. ^ model 3 is the case in our MSSM sam- 
ple which produces the flattest spectrum and therefore 
may induce the most severe break. It is a very heavy neu- 
tralino, = 2730 GeV, and the rescahng needed in this 
case may be uncomfortably high, /(5 ~ 10^ (but still it 
does not violate any experimental constraint from other 
cosmic ray measurements). The second example, model 4, 
has a mass of 1400 GeV and f5 = 7200. Clearly these 
are extreme cases: in this category we may include a large 
fraction of models in our sample, which with a very mild 
rescaling give a distorsion of the flux at the highest en- 
ergies. The chance of detecting less severe breaks in the 
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Fig. 5. Estimate of possible fluxes in case of finite an- 
tiproton lifetime Tp, as compared with the standard back- 
ground case when Tp = oo. The total flux is the sum of a 
neutralino induced component (not shown) and the sec- 
ondary component (dashed lines) which, as plotted in the 
3 examples, can be much suppressed compared to line b. 

energy spectrum clearly depends on how accurately the 
antiproton flux is measured. 

3.3. Finite antiproton lifetime 

In exotic scenarios, the antiproton may have a finite life- 
time Tp. In the cosmic ray context this is an interesting 
ingredient in case Tp is lower than the characteristic escape 
time of antiprotons from the Galaxy. The escape time is 
determined by the parameters which define the propaga- 
tion model; in our case it is about 10 Myr. This value is 
clearly much lower than the bound one can infer in CPT 
conserving theories from the experimental limit on the 
proton lifetime {Tp = Tp > 10'^^ yr, see Caso et al. (1998)). 
On the other hand it is much larger than the most strin- 
gent direct experimental bound (r^ > 0.05 Myr from the 
p — > fi~j decay mode, see Hu et al. (1998)). 

Geer & Kennedy (1998) have recently claimed that 
from cosmic ray measurement it is possible to infer a limit 
of Tp > 0.8 Myr (90% C.L.). In their analysis the sec- 
ondary component only is included. Adding a neutralino 
signal obviously may change this result. Actually we show 
here that in a clumpy scenario, even in case of a finite an- 
tiproton lifetime, we can produce models which give at the 
same time a very good fit of the existing Bess data and 
whose spectral features are peculiar enough to be distin- 
guished from the standard background once measurements 
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at higher energies will be available. In Fig. ^ we compare 
three such models with the case of standard background 
and infinite Tp. Model 5 is a 51 GeV gaugino-like neu- 
tralino whose antiproton flux has been scaled by / <5 = 49 
and for Tp = 0.82 Myr. Also shown in the figure is the 
reduction in the background flux induced by considering 
Tp = 0.82 Myr (dashed line labelled by 56). Model 7 is 
the heaviest model for which the predicted flux is still 
in excellent agreement with existing data (90% C.L. fit 
for = 477 GeV, f S = 1.2 ■ lO'^ and Tp = 2.92 Myr), 
while model 6 is some intermediate case (m^ = 188 GeV, 
f S ^ 78 and Tp — 1.32 Myr). The trend is that for heavier 
neutralinos there is a larger overproduction of antiprotons 
in the high energy range. 

Applying the largest possible rescalings consistent with 
7-ray measurements we find that a 56 GeV neutralino 
model gives a flux which is consistent with Bess data 
at 90% C.L. in case the antiproton lifetime is as low as 
Tp = 0.15 Myr. The bound of Geer & Kennedy (1998) is 
clearly violated. Notice that we are comparing with a more 
aboundant data set than in that reference (Bess 97 data 
were not included there) and that we used our standard 
values for the parameters which define the diffusion model 
and solar modulation. If uncertainties were included the 
lower bound we would get with this method would proba- 
bly be very close to the most stringent direct experimental 
bound Tp > 0.05 Myr or lower. 
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4. Conclusion 

To conclude, we have shown that there is a chance of de- 
tecting neutralino dark matter in upcoming measurements 
of the cosmic antiproton flux at high energies. The signa- 
tures we propose here are alternative to the signature of 
an exotic component at low kinetic energies, which seems 
not to be required by present data. We have also discussed 
the possibility that antiprotons have a finite lifetime and 
shown that the limit on Tp which is possible to set on the 
basis of cosmic ray measurements is comparable to those 
in direct experiments. 
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